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Abstract. We consider a zero-surface-tension two-dimensional Helc-Shaw flow in an infinite wedge. There 
exists a self-similar interface evolution in this wedge, an analogue of the famous Saffman- Taylor finger in 
a channel, exact shape of which has been given by Kadanoff. One of the main features of this evolution 
is its infinite time of existence and stability for the Hadamard ill-posed problem. We derive several exact 
solutions existing infinitely by generalizing and perturbing the one by Kadanoff. 



1. Introduction 

The Helc-Shaw problem involves two inmiscible Newtonian fluids that interact in a narrow gap between 
two parallel plates. One of them is of higher viscosity and the other is effectively inviscid. The model under 
consideration is valid when surface-tension effects in the plane of the cell are negligible. In the most of 
the cases it is known that when a fluid region is contracting, a finite time blow-up can occur, in which a 
cusp forms in the free surface. The solution does not exist beyond the time of blow-up. However, Saffman 
and Taylor in 1958 |13j discovered the long time existence of a continuum set of long bubbles within a 
receding fluid between two parallel walls in a Hele-Shaw cell that further have been called the Saffman- 
Taylor fingers. It is worthy to mention that the first non-trivial explicit solution in the circular geometry 
has been given by Polubarinova-Kochina and Galin in 1945 I12j . They also have proposed a complex 
variable approach, that nowadays is one of the principle tools to treat the Hele-Shaw problem in the plane 
geometry (see, e.g., [8 [ 114) 1 . Following these first steps several other non-trivial exact solutions have been 
obtained (see, e.g., [2l l3l HI 151 l9l llOllllj ). Through the similarity in the governing equations (Hele-Shaw 
and Darcy) , these solutions can be used to study the models of saturated flows in porous media. Another 
typical scenario is given by Witten-Sandcr's diffusion- limited-aggregation (DLA) model (see, e.g., [T]l. 
In both cases the motion takes place in a Laplacian field (pressure for viscous fluid and random walker's 
probability of visit for DLA). One of the ways, in which several new exact solution have been obtained, is to 
perturb known solutions. For example, Howison [7j suggested perturbations of the Saffman- Taylor fingers 
that led him to new fingering solutions keeping the same asymptotic behavior as time t — > oo. Recently, 
Hele-Shaw flows and Saffman- Taylor fingering phenomenon have been studied intensively in wedges (see, 
e.g., IU [2j G3 El HOI IH| nad the references therein). In particular, Kadanoff |10| suggested a self-similar 
interface evolution between two walls in a Helc-Shaw cell expressed explicitly by a rather simple parametric 
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function with a logarithmic singularity at one of the walls. By this note we perturb Kadanoff's solution 
and give new explicit solutions with similar asymptotics. 

2. Mathematical model 

Wc suppose that the viscous fluid occupies a simply connected domain fl(t) in the phase z-plane whose 
boundary T(t) consists of two walls and r 2 (£) of the corner and a free interface T 3 (t) between them 

at a moment t. The inviscid fluid (or air) fills the complement to fl(t). The simplifying assumption of 
constant pressure at the interface between the fluids means that we omit the effect of surface tension. The 
velocity must be bounded close to the contact point that yields the contact angle between the walls of 
the wedge and the moving interface to be ir/2 (sec Figure A limiting case corresponds to one finite 
contact point and the other tends to infinity. By a shift we can place the point of the intersection of the 
wall extensions at the origin. To simplify matter, we set the corner of angle a between the walls so that 
the positive real axis x contains one of the walls and fix this angle as a £ (0, ir]. In the zero-surface-tension 



iy 




Figure 1. il(t) is the phase domain within an infinite corner and the homogeneous 
sink/source at oo 



model neglecting gravity, the unique acting force is pressure p(z, i) = p(x, y, t). The velocity field averaged 
across the gap is given by the Hele-Shaw law (Darcy's law in the multidimensional case) as V = — Vp. 
Incomprcssibility implies that p(z,t) is simply 



(1) 

The dynamic condition 
(2) 



Ap = 0, in the flow region fl(t). 



P 



0, 



is imposed on the free boundary T3 = T3 (t) . The kinematic condition implies that the normal velocity v n 
of the free boundary T3 outwards from Q (t) is given as 

dp 

on r 3 



(3) 
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On the walls I\ = Ti(t) and T 2 = ^(i) the boundary conditions are given as 

(impermeability condition). We suppose that the motion is driven by a homogeneous source/sink at infinity. 
Since the angle between the walls at infinity is also a, the pressure behaves about infinity as 

p ~ — — log \ z\, as I z I — > oo , 
a 

where Q corresponds to the constant strength of the source (Q < 0) or sink (Q > 0). Finally, we assume 
that ^(0) is a given analytic curve. 

We introduce the complex velocity (complex analytic potential) W(z,t) = p(z,t) + iip(z,t), where 
— tp is the stream function. Then, Vp = dW/dz by the Cauchy-Ricmann conditions. Let us consider 
an auxiliary parametric complex £-plane, ( = £ + irj. We set D = {£ : |£| > 1, < arg<^ < a}, 
D 3 = {z : z = e ie , 6 G (0,a)}, D 1 = {z : z = re la , r > 1}, D 2 = {z : z = r, r > 1}, 91? = D 1 U D 2 U L> 3 , 
and construct a conformal univalent time-dependent map z = /(£, t), f : D — *• O(t), so that being continued 
onto 3D, /(oo, t) = oo, and the circular arc -D3 of dD is mapped onto T 3 (see Figure|2J). This map has the 




Figure 2. The parametric domain D 



expansion 

00 

/(C,t) = (E a ^)C" ? 

about infinity and Oo(t) > 0. The function / parameterizes the boundary of the domain Q(t) by Tj = {z : 

ai = /(C,t),Cei? j },i = l,2,3. 

We will use the notations / = df/dt, /' = df/d£. The normal unit vector in the outward direction is 
given by 

/' „ . i 

h = — CtjTT 011 1^3, n = —i on T2, and h = ie la on Y\. 
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Therefore, the normal velocity is obtained as 
(5) v n — V • h - 



- Re {^W\ ] 



dn I 0, forCG^i 
0, for (eD 2 

The superposition W o f is the solution to the mixed boundary problem QJ, ©, in D, therefore, it is 
the Robin function given by Wo f = — 2 log£. On the other hand, 

r Re (/Cf/l/'l), forced 

(6) i>„ = < Im (/e~ M ), for (eA 

1 -Im(./), forCG-D 2 

The first lines of ® give us that 

(7) Re(fCf) = -, forCe^. 

a 

The resting lines of J^J, © imply 

(8) lm(/e- m )=0 forCe£»i, Im (/) = for £ G -D 2 - 

3. Exact solutions in a wedge of arbitrary angle 

We are looking for a solution in the form 



\— C + C3(C), 



where 5(C) is regular in D with the expansion 

00 

.9(0 = £ 5r 

about infinity. The branch is chosen so that <?, being continued symmetrically into the reflection of D is 
real at real points. The equation Q implies that on D3 the function g satisfies the equation 

Re (g(C) + C<?'(0) = 0, C e D 3 . 

Taking into account the expansion of g we are looking for a solution satisfying the equation 
(9) 0(C) + Cs'(C) = p^, (£D. 

Changing the right-hand side of the above equation one would obtain other solutions. The general solution 
to © can be given in terms of the Gauss hypergeometric function F = as 

Cfl(C) = C-2CP(-,i,i + -;-C 5 )+c. 



We note that /' vanishes for = (2/(1 + y/2Qt/a)) — 1, therefore, the function / is locally univalent, 
the cusp problem is degenerating and appears only at the initial time t = and the solution exists during 
infinite time. The resulting function is homeomorphic on the boundary dD, hence it is univalent in D. 
This presents a case (apart from the trivial one) of the long existence of the solution in the problem with 
suction (ill-posed problem). To complete our solution we need to determine the constant C . First of all we 
choose the branch of the function 2F1 so that the points of the ray £ > 1 have real images. This implies 
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ht ( a ) ( b ) 

FIGURE 3. Interface evolution in the wedge of angle: (a) a = 27r/3; (b) a = it/3 

that Im C = 0. Wc continue verifying the asymptotic properties of the function f(e ,t) as 6 — > a — 0. 
The slope is 

lim a,Tg[ie l6 f'(e ie ,t)} = a + £ + lira arg [ J — + ^ ) = a + tt. 

e^a-o 2 e^a-o \ V a e « + 1 / 

To calculate shift we choose C such that 

lim Im [e- ia f'(e ie ,t)] = 0. 

Using the properties of hypergeometric functions we have 

lim ImF(-,l,l + -;e i7 l = -. 

7^0+0 \7T 7T / 2 

Therefore. C — a. We present numerical simulation in Figure [3] 
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Figure 4. Kadanoff's solution 
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The special case of angle a = ir/2 has been considered by Kadanoff |10| . The hypergeometric function 
is reduced to arctangent and we obtain 

1 +l( 7T 



(10) 



/(C,*) = (V4QiA + i)C + *iog 



Q > 0. 



1 - i( 2 

This function maps the domain > 1, < arg ( < tt/2} onto an infinite domain bounded by the 
imaginary axis (ri), the ray T 2 = {r : r > ^AQt/n + 1} of the real axis and an analytic curve which 
is the image of the circular arc, see Figure 0] 

4. Perturbations of Kadanoff 's solution 

Kadanoff 's solution (jlC)|t can be thought of as a logarithmic perturbation of a circular evolution with 
the trivial solution /o(C *) = V 4Qi/7r£. A simple way to generalize the solution Hl()|l is to perturb another 
function. For example, one may choose 



oi, A = J^i 



Wc find the solution f((,t) in the form /(C,£) = fo((,t) + ^(C) similarly to the preceding section. Then 
the equation is satisfied when 



Re 



C*'(C) 



0, or Re 



= 0, 



/<>(£*) cC + l/cC 

where h! ~ (£ — i) as £ — > i in the unit circumference. Wc choose a consistent form of /i as 

cC 2 V(C) _ c 2 - 1 



Integration yields 



c 2 C 2 + 1 C 2 



fc( o = cc + l-i( c -i)iog^±i 



c, 



where C is a constant of integration. Satisfying the conditions on the walls we deduce that C = 0, and 
finally, we get a logarithmic perturbation of the elliptic evolution as 



f(C,t) = (AVt + l) (<+£) -<(c- 
see the interface evolution in Figure 03 



log 




Figure 5. Logarithmic perturbation of the elliptic evolution 
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The next goal is to obtain perturbations of the logarithmic term of Kadanoff 's solution that with the 
same asymptotic as t — ► oo, such that the interface has finite contact points at a finite moment. Let us 
consider the function 

ff(C,t) = 2d(t)C-log^Ct) 
C-a(i) 

The functions a(t), d(t) are to be chosen such that equation (JZJ) is satisfied for the moving interface as 
well as the conditions of impermeability and univalence hold. The local univalcncc is followed from the 
first restriction a(t)/d(t) < 1. Substituting H into equation Q and comparing the Fourier coefficients, we 
derive the following system of equations for the functions a(t) and d(t): 

(1 + a 4 )dd - a 3 d + da 2 a -aa = ^L(l + a 4 ), 
-2a 2 dd + ad- ad = -3-a 2 

This system can be easily solved and the first integrals are 

(ii) d(t) = 



1+ v/l + 4a 2 (*)(QtA-Ci) 



2a(t) 



(12) 
where 



m l + a 2 (t) 
a(t) 8 l-a 2 (t) 



C 2 



o(0)' ~a(0) '"°l-a 2 (0)' 

are the constants of integration. Let us assume the initial condition a(0) € (0, 1). Making use of the system 
(|11I12J) we arrive at the explicit function t(a) inverse to a(t) 



(13) 



*(-) - Q 



a 2l 0g l±^ 



+ a 2 C 2 



4a 2 



that exists, is continuous, and increases in the interval a € [a(0), 1). Therefore, the function a(t) increases 
from a(0) to 1 as t £ [0, oo). By ltTT|l we conclude that d(t) ~ 0{\fi) as i — » oo. The rotation of is 
exactly Kadanoff 's solution when a = 1, and d(t) is appropriately chosen as in l|10(l . To make a numerical 
simulation one may use the Newton method of the solution of a non-linear system (Howison [7] presented 
the numerical approximation of an analogous solution in a narrow channel), see Figure Choosing a(0) 




Figure 6. Perturbation of Kadanoff 's solution for a(Q) = 0.5, d(0) = 1 
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rather close to 1, one may give an explicit analytic approximation by, e.g., introducing two functions 



d(t) = 



l+y/l + 4(Qt/n-Ci) 



f exp(2rift) - C 2 ) - 1 
exp(2d(i) -C 2 ) + l' 



The initial conditions a(0) and d(0) are to satisfy the inequalities 1 — 4Ci > 0, 2d(0) — C2 > 0. To proceed, 
we simplify putting a(0) = d(0). Then these inequalities are satisfied for d(0) € (^/3/4, 1). It is easily seen 



from lO that |1 - a{t) \ ~ e - A Then \d(t) - d(t)\ ~ too. Similarly, ~ y/t and |1 - a(i)| 
Both a(i) and a(i) tend to 1 rapidly and the error \d(t) — a(t)\ is of the same order for t ~ 00. 
Now we evaluate the error |d(t) — a(t)\ for < t < 00, and claim that 

(14) 0<a(£)-o(£)<8(a(0)-6(0)). 

To prove this we estimate the distance p(a) between the inverse function (|13|l and 

/ 



t(o) 



7T 



C 2 



V 



as 



p(a) = i(a) - t(a) = (1 
The derivative of p is 

p'(a) = 2a (log 



log' 



C, 



Ci 



and /9«0)) = 2(1 - a(0)). 



1 - a 2 



Co 



log 



1-a 2 



Since the function log increases, the function p(a) may have a critical point a c , a c G [a(0), 1), which 

is the maximal solution to the equation 1 _ ) 4 a2 = log yz^i — C2 in the interval [a(0), 1). The latter equation 
implies 



p(a c ) = (1 - a 2 ) log 



C 2 = 



16(1 -al) 
(l + « 2 ) 2 



<16(l-a(0) 2 ) = 8p(a(0)), 



that proves (|14f) . 

Moreover, a(0) — d(0) is decreasing and non-negative as a function of the initial condition a(0), that 
vanishes as a(0) — * 1. Therefore, given a small positive number e, we may choose a(0) = d(0) close to 1 
such that d(t) approximates a(t) with the precision e during the whole time < t < 00. Desired quantity 
a(0) satisfies the equation 



a(0) 



exp( v /4 a (0) 2 -7 + log|±gf -l) 



exp 



4a(0) 5 



7 + log 1+a W 2 

' ' 10 8 l-a(0) 2 



A similar conclusion may be made for the function d(t) and its approximation d(t) (note that d(t) < d(t) < 
1). Moreover, the mapping 

c + m 



H(C,t)=2d(t)C-log 



C-a(i) 



converges to Kadanoff 's solution as t 
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FIGURE 7. Error \a(t) - a(t)\, for a(0) = 0.9, d(0) = 0.9 
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